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f/3 I A theoretical derivation is given for the formula describing A^-electron ion- 

ization of atom by a dc field and laser radiation in tunneling regime. Numerical 
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J>-,. examples are presented for noble gases atoms. 
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1 Introduction 

Many-electron ionization of atoms by laser field was firs observed by Suran and Zape- 
sochny ^ in alkali-earth atoms (the review of that work as well as some earlier ones 
see in 10). At present, such studies form one of main guidelines in physics of strong 
field interaction with atoms [^. 

A number of theoretical models were proposed for interpretation of the gathered 
experimental data. Some models dealt with direct influence of laser radiation on the 
atomic electrons [^ ||, ^ |^, ||, the others consider highly- stripped ion formation due 
to nonelastic scattering of previously emitted electrons with the parent ion [^ ^, [11 
These models allow to explain a number of observed features of the phenomenon ||12|, 



|T^, [T^, . Nevertheless, there are some difficulties in theoretical description of highly- 
stripped ion formation in laser field which is not related to nonelastic collisions |Q, |^, 
1^, 0, ||. Due to these difficulties, the above mechanisms cannot be properly used for 
explanation of the experiment. 

At the same time, it is well-known fact that the single-charged ion formation by 
a laser field in tunnelling regime can be satisfactory described in terms of relatively 
simple formulae of the ADK theory [ll7| , |18|, 0. An empirical generalization of the 
ADK formulae for describing the highly-stripped ion formation was proposed in |2^ . 



So it would be reasonable to generalize the available theory of tunnelling in atoms to 
the case of non-sequentional multiple ionization of atom. Solution of this problem is 
the objective of the present work. 

Obviously, the Josephson effect can be considered as a solid-state analogue of the 
considered phenomenon. Some considerations on difference between the one- and 
many-particle tunnelling are mentioned in reference |^. Comparison of these con- 



siderations with the results of the present work shows that the mentioned difference 



for tunnelling in atoms is not so trivial as it was described in [Bll . 



2 Asymptotics of the many-electron wave function 

Let us remembers some facts which make the main proposed concepts easer to un- 



derstand. To describe optical transitions in complex atoms. Bates and Damgaard [22 



modified the Slater method p3|. Basically, the nodeless character of Slater orbitals 



was retained. Unlike the Slater method, the effective nuclear charge ceases to be a 
fitting parameter for valence electrons in atom since it coincide with the residual ion 
charge. But the effective principal quantum number is uniquely determined by the 
electron coupling energy. So the asymptotical region of electron motion is considered, 
where the atomic potential has Coulomb shape. High accuracy of oscillator strengths 



calculations P3 using the Bates-Damgaard method and its clear physical justification 



allows to use this method in calculations of other atomic characteristics determined by 
large electron-nucleus distances. 

The tunnelling probability is also determined by large electron-nucleus distances 
where the energy of the electron interaction with the external field becomes compa- 
rable with the attractive energy of the residual ion. So the Bates-Damgaard method 
can be used for describing the tunnelling effect. Such a procedure was developed in 
recent work p5[ for tunnelling calculation in Rydberg molecules. In that work some 
evaluations are presented for the applicability conditions of the method. 

Let A^ equivalent (i. e. belonging to the same atomic shell) electrons are removed 
from the atom via tunnelling. Then the asymptotic behaviour of the radial part of 
A^-electron wavefunction in the Bates-Damgaard approximation is determined by the 
product of properly symmetrized one-electron function asymptotics: 
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C.i = (27rz.)-V2 (^IJ Lie), L{e) = ^^^j (1 - e^)-/^. 

Here b = au/Z, Z is the residual ion charge, a = h'^/fie'^ is Bohr radius, /i, e are the 
mass of electron and the absolute value of its charge, e = {I + \)/v- The C^i constant 
in (|l|) is determined in quasiclassical approximation not implying the condition / -C v., 
which was required in |19[. It results in the arising of L{e) function with L{e) — * 1 
at e — i> 0. After passage to this limit the expression (|l]) for the C^i constant turns 
into the formula (11) of the reference |l^ (with an inaccuracy corrected: the number 
e = 2.718 . . . should be omitted). 

The expression (|l]) for C^i is obtained under e < 1. For e > 1, the quasiclassical 
approximation is not valid, so calculation of C^i requires numerical approaches (see, 
e. g. 
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The principal quantum number u is determined by the electron coupling energy. 
Denoting the first, second etc. ionization potentials of the atom as -Ei/e, E2/e . . ., the 
principal quantum number of j-th removed electron is 






' \Z^e\ 

If the electron are equivalent and are simultaneously removed from the atom, then for 
all the electrons 

(2) 



NZ^e^ 
where 

AT 

En = ^ Ej 



is the coupling energy of A^ electrons. Note that in framework of the considered model, 
the asymptotic behaviour of the bound electron wave function (|I|) depends on the 
number of the removed electrons. So a partial account is provided for many-electron 
effects in the initial state. 

Now we consider A^-electron ionization as removal of a A^-electron "bundle" - a 
peculiar kind of quasiparticle of mass Nfi and of charge —Ne. In the region which 
determines the ionization process, we consider the distances between the electrons in 
the bundle to be much less than the separation between the atomic core and the center 
of bundle mass. Denoting the distance between the i-th and j-th electrons as Xij, and 
the position of the center of bundle mass as R, we write the corresponding inequality: 

Xij < R. (3) 

Since the atom-laser radiation interaction is considered in dipole approximation, the 
influence of the field on A^ individual electrons is completely equivalent to the influence 
of the field on a quasiparticle of charge — A^e which is located at the point R. As for 
the interaction of this quasiparticle with the core Coulomb field, the correspondent 
error value is ~ {xij/RY, which is small due to the accepted inequality @. 

For the mathematical description of the considered model, one should solve a prob- 
lem which is analogous to that is occurred, e. g. in nuclear a-decay theory. This 
problem is to construct the quasiparticle wave function '^\^,{R, {xi}) at large dis- 
tances from the residual system, using the one-particle wave functions of the system 
in the initial state. Symbols in the braces are sets of quantum numbers or coordinates 
of individual particles. To solve this problem we consider the asymptotics of the func- 
tion ^E'L/mj at -R ^ oo, which is a product of the one-electron function asymptotics (|I]). 
It is easy to see that the radial dependencies of the functions (|1|) bring the factor 

NR\ fRY'^""'^ 



into the asymptotics of ^L/^}- '^^ obtain the angular dependence, the mean of the 
R, {xi} variables should be detalized. Since the problem has the axial symmetry for 
the linearly polarized field, the orbital moment projections of non-interacting electrons 
onto the polarization direction are conserved. So it is convenient to leave the azimuth 
angles (pi the same that in the original spheric coordinate system centered in the atomic 
nucleus. The change of variables will effect only on the absolute values {tj} and polar 
angles {Oi}. At 6 ^ 0, the behaviour of the Legendre polynomials involved in the 
spheric functions (|ID, is determined by 

P^I(cos^,) ~ (-1)H!^L^ = (_i)Hi!i il^l-- 



Substituting here ri -^ R, Vi^ -^ Rz and introducing the parabohc coordinates ^ = 
R + Rz, 1] = R — Rz for the center of the bundle mass, the asymptotics of the A^- 
electron function at ^ ^ rj can be written in the form 
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(f){^, 7]) ~ exp 
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(4) 



Here x is the normahzed per unit wave function of the electron inner motion in the 
bundle. Note that there are only 2{N — 1) independent variables {ri,6i} of 2N. The 
function (/)(^, rj) describes the motion of the center of the bundle mass. 



3 Tunnelling probability 

The further calculation of the tunnelling probability is implemented according the 
standard technique [1^, |2^, an account provided for that the electron bundle mass 
is Nfi and its charge is —Ne. Substituting the function 0(^, 77) from (^ into the 
Schrodinger equation 

d fA4>\ ( ^ E^Nii 



e- 
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0, 



de V d^y V 2/12 

describing the motion with respect to the parabolic ^ coordinate at ,^ ^ 00, we obtain 
the variables separation constant: 

M -I 



^^"t 



iV(z/-l) 



(5) 



The centrifugal potential is neglected since it vanishes rapidly at ^ — i> cxd. 

Now we consider the external field F{t) to be slow-varying, and use quasiclassical 
approximation for the wave function (f)F{^,v) which describes the center of the bundle 
mass motion in the field. In the below-threshold domain 

1 r( 
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where ^1 is the greater root of the equation p(^) = 0. Comparing the expression 
with the function 0(^, 77) from (Q) at the point ^0 lying in the region 
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we obtain the x value: 

1/2 / T ri 



x{mio) ^ (^) exp (-ly W)\ddj <P{^o,v)- (9) 

The condition of existence of the region (P) leads to following restriction to the external 



field: 



eZ e ( Z^ 



F«F, = — = - - , (10) 



which differs from the condition arising in the one-electron tunnelling description only 
by the definition of the v value. It should be noted that for v essentially greater than 



1 (what holds, e. g. for Rydberg states) the inequality ([T0|) is changed by a stronger 
one: 



which is deduced from the condition of existence of the potential barrier [^ . 

The formulae (^) and (P) determine the function (J)f{^,v) outside the barrier. With 
the account of inequality (H), its squared absolute value is ||2j 

urn, /e.y'"-"/.,^-" 
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Using (0) and (^), it is easy to see that the dependence on the arbitrary parameter ^q 
is actually disappeared in ([T^): 



The ionization probability is determined by the flux of probability density ( |T3l) through 
a plane perpendicular to 2;-axis ||27[] : 



H^if)(F)~27r^~t;,|0^(e,r^)|Vdp, ^^ = ^' P = V^, dp ^ J^dr/. 



Substituting here the formulae (H) and (|T3|), we obtain: 



TThM\{2l + l)''C^^ fzV^-^^ il + \mj\)\ 



"^ ^ ' a?^ji 2*^-2ArM+i y^, ) 11 (|^^.|!)2(/ _ |^^.|)! 

2^^^2A.(._i)-Am . 2iVF,\ 

exp (^-^^ j . (14) 

This formula determines the A^-electron tunnelling probability in dc field within a factor 
accounting for the overlapping of wave functions of the electrons remaining in the atom, 
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with the wavefunctions of the same electrons in the initial state. Obviously, this factor 
cannot exceed 1, and its more accurate evaluation can be performed only numerically. 
Note that the A^ multiplier in the exponent in (jT^) in now ways gives an exhaustive 
account for the dependence of this exponent on A^, as it was considered in |^. Due 
to the formulae (0), (^^, this dependence is significantly more complicated and it is 
determined by the spectrum of the particular atom. We present below (figure |1|) a 
numerical example illustrating this statement. 




Figure 1: Relation of the 3-charged ion formation probabilities for noble gases atoms 
by two different channels (see the text). 

Now we consider that 

F{t) = Fo cos iut, (15) 

where u is the laser field frequency. It is a well-known fact that the tunnelling in a 



laser field is possible for small values of the Keldysh parameter |29 
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where Ei is the coupling energy of one electron. Following the technique developed 
2^] for "particle" of mass Nfi and charge —Ne, it is easy to see that the iV-electron 



m 



tunnelling is possible for small values of the parameter 



In 



eF 



-u. 



(16) 



Since the coupling energy is increasing for each subsequent electron, A^-electron tun- 
nelling requires field values lower than A^-electron tunnelling cascade. 

Substituting (|15|) into (0), we average the result over the time interval t G 
[— 7r/2to', 7r/2u;] []18|[]. Due to the inequality (p!0[), the integral arising here can be cal- 
culated using the saddle-point method. Under the condition (ITT]) fulfilled, the saddle 



^The values t E [Tr/2uj, 3tt/2lu] leads to F{t) < and the tunnelling takes place in the direction of 
negative z semiaxis. 
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point is t = 0, and the final formula is: 



^hM\{2l + l)^C^f fzV^-^^ il + \mj\)\ 



"^ '^ O; 2,, 9Af-3/2 ArM+3/2 \ ,, 11 



a2/i 2A'^-3/2iVM+3/2 \^^^y ii(|^^.|!)2(/_|^^.|)! 



Remember that the exponent dependence on N in (|T^) is not reduced to the factor A^ 
which is written explicitly. 

4 Numerical examples 

Unfortunately, the obtained formulae cannot be immediately related to an experiment, 
because, along with the direct A^-fold ions formation, there are a number of cascade 
processes as well as other ionization mechanisms due to nonelastic collisions of electrons 
and ions p, |10|, |lT|. For the relation of the theory with an experiment, the correspon- 



dent kinetic equations are to be solved, that should be a subject for another work. So 
only some illustrative examples are considered in this section. 

The figure |l] presents the relation of probabilities of 3-fold ions formation in 
the noble gases resulted from two 2-cascade processes: A -^ A+ -^ A'^"*" and 
A -^ A2+ -^ A^^. These probabilities are denoted as W{1;2) and W{2;1) corre- 
spondingly. They have similar dependence on the laser pulse duration. As it is seen, 
the relation W{1] 2)/W{2; 1) is not equal to 1, as it is follows from the results of refer- 
ence 
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The following result seems to be curious. The 2-electron tunnelling probabili- 
ties for neutral atoms can be greater than the one-electron tunnelling probabilities 
in correspondent singly charged ions. E. g., for Ar atom the 2-electron tunnelling 
probability exceeds the 1-electron process probability for Ar"*" ion at the intensi- 
ties I > 10^^'^^W/cm . The same result takes place for Kr at / > 10^^'^^W/cm , 
for Xe at / > 10^^'^^ W/cm . At the same time, for light noble gases atoms He and Ne, 
the probabilities of one-electron tunnelling in singly charged ions are approximately 
by two orders greater than the probabilities of two-electron process in the correspon- 
dent neutral atoms at / ~ 10^^ W/cm . These facts shows wide range of experimental 
situations arising in multiphoton tunnelling effect. 
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